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Abstract: First-principles calculations based on density functional theory and the pseudopoten- 
tial method have been used to investigate the energetics of H2O adsorption on the (110) surface 
of TiC>2 and SnCV Full relaxation of all atomic positions is performed on slab systems with peri- 
odic boundary conditions, and the cases of full and half coverage are studied. Both molecular and 
dissociative (H2O — > OH~ + H + ) adsorption are treated, and allowance is made for relaxation of 
the adsorbed species to unsymmetrical configurations. It is found that for both Ti02 and SnC-2 
an unsymmetrical dissociated configuration is the most stable. The symmetrical molecularly ad- 
sorbed configuration is unstable with respect to lowering of symmetry, and is separated from the 
fully dissociated configuration by at most a very small energy barrier. The calculated dissociative 
adsorption energies for TiC>2 and SnC>2 are in reasonable agreement with the results of thermal 
desorption experiments. Calculated total and local electronic densities of states for dissociatively 
and molecularly adsorbed configurations are presented and their relation with experimental UPS 
spectra is discussed. 



1 Introduction 



Most oxide surfaces react readily with water and become partially covered with molecular H2O 
or hydroxyl groups on exposure to air under ambient conditions. This adsorption of water has 
important consequences for surface processes such as catalysis and gas sensing. It is widely 
accepted that molecular water will generally bind to most oxide surfaces through the attraction 
of its electric dipole to the ionic charges. However, the mechanisms by which dissociative ad- 
sorption occurs are still very poorly understood. It has sometimes been suggested that surface 
defects such as vacancies or steps are essential to water dissociation, and there is no doubt that 
there are some oxide surfaces for which this is true. For example, recent calculations using both 
semi-empirical [0 and first-principles techniques have shown that water dissociation is not 
energetically favorable on the perfect MgO (100) surface, but that it is favorable at steps and 
corners. However, it is unlikely that defects are essential for dissociative adsorption in general, 
since oxide surfaces vary enormously in their geometry and electronic structure. 
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The aim of this paper is to present first principles calculations on the energetics of H2O adsorp- 
tion on the SnC>2 and Ti02 (100) surfaces. Both materials have the rutile crystal structure, and 
the (110) surface is the most stable when the materials are stoichiometric. The stoichiometric 
(110) surface is far from flat, so we are dealing with a situation which is geometrically very dif- 
ferent from the MgO (100) surface. Although they have the same geometry, SnC>2 and TiC>2 are 
electronically different, since Ti is a transition element and Sn is a group IV element. Our hope 
is therefore that by comparing with previous work on the interaction of H2O with MgO (100) 
we can learn about the influence both of geometry and of electronic structure. Our choice of 
Sn0 2 and Ti0 2 is also motivated by their technological importance. Sn0 2 is widely used in 
gas-sensing devices, and adsorbed water is known to have an important effect on the adsorp- 
tion of other molecules ||]. The interaction of Ti02 with water is important in photocatalytic 
applications 0. 

A considerable experimental effort has gone into the study of H 2 adsorption on Ti02 (HO), 
though a clear picture has yet to emerge. Hydroxyl groups on the surface have been detected 
at 300 K in ultraviolet photoelectron spectroscopy (UPS) experiments by Henrich et al. ||, 
and their presence was later confirmed in synchrotron radiation studies by Madey's group |7|||. 
The amount of dissociated water was less than a monolayer, and was weakly dependent on 
the oxygen vacancy concentration before adsorption. The minor influence of point defects and 
steps on the reactivity with water has also been reported by Muryn et al. |J. Experiments us- 
ing a combination of thermal desorption spectroscopy (TDS), x-ray photoelectron spectroscopy 
(XPS) and work-function measurements have recently been reported ]l(J. A three-peak desorp- 
tion spectrum was found, with features at 170 and ~275 K and a high-temperature tail up to 
~375 K. The two low-temperature features appear to be due to molecular water and the high- 
temperature tail to the reaction of surface OH groups to form H 2 0. A widely accepted model 
is that of Kurtz et al. J7|, who proposed that adsorption on nearly perfect Ti02 (HO) occurs 
by molecular adsorption of H 2 at 5-fold coordinated cation sites, followed by dissociation to 
give OH~ attached by its oxygen end to the cation, and H + bonded to lattice oxygen to form a 
second type of hydroxyl group. The presence of at least two types of adsorbed OH~ has been 



indicated by infrared spectroscopy [11,12 



The experimental situation for Sn02 is rather similar, although the experimental evidence 



is less extensive. Recent experiments by Gercher and Cox (13[ using a combination of TDS 
and UPS show the presence of disso datively adsorbed water on Sn0 2 (HO). On the perfect 
stoichiometric surface, three distinct TDS peaks are seen at T = 200, 300 and 435 K. The 200 K 
peak is attributed entirely and the 300 K peak mainly to molecular H 2 0. The feature at 435 K 
is assigned to disproportionation of surface OH~ to form H2O. As in the case of Ti02, surface 
defects do not appear to play a major role. 

Few theoretical investigations of H2O adsorption on rutile surfaces have been reported. Jay- 



cock and Waldsax [14] performed calculations for the Ti02 (HO) and (100) surfaces using an 
empirical interaction model, which predicted adsorption energies much greater than the ex- 
perimental values. The dissociative adsorption of single molecules on clusters embedded in a 
point-charge lattice has been studied using the Xa approach for Ti0 2 (HO) by Tsukada et 
al. WM. The dissociative adsorption of single water molecules on Ti02 as well as the case of full 
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hydroxylation have been studied using a semi-empirical technique by Goniakowski et at [|T(J 
and Goniakowski and Noguera [[jj]. Both molecular and dissociative adsorption of water on 
Ti02 have been treated using a periodic Hartree-Fock approach by Fahmi and Minot ||17|| , who 



find a strong preference of water to dissociate. The above studies considered only very sym- 
metrical configurations of adsorbed species, and all found that both molecular and dissociative 
adsorption are energetically favorable. The only attempt to study relaxation of surface species 
to unsymmetrical configurations has been the very recent semi-empirical study by Bredow and 
Jug on the adsorption of H 2 on rutile and anatase Ti0 2 surfaces [Q. To our knowledge no 
theoretical attempts to model water adsorption on Sn0 2 surfaces have been reported. 

The present work is based on density functional theory (DFT) in the pseudopotential ap- 
proximation. This approach has been widely used for oxides, including MgO fL9| , Li 2 |2(J 



Al 2 03 PH , Ti0 2 p2| , p3| and Sn0 2 p4| , and is known to give accurate results for the energetics 



of perfect crystals, lattice defects, surfaces and molecular adsorption. We have recently reported 
a detailed study of the stoichiometric and reduced Sn0 2 (110) surface using this approach f25fl . 
However, since the energetics of dissociation is not generally described very accurately within 
the standard local density approximation (LDA), we include gradient corrections in the present 
work, using the Becke-Perdew scheme |2B],|27| . We present here results on the energetics and re- 
laxed geometry of both molecularly and dissociatively adsorbed water on Sn0 2 and Ti0 2 (110), 
and we study both symmetrical and unsymmetrical modes of adsorption. In addition, we report 
total and local electronic densities of states, through which we attempt to make contact with 
recent spectroscopic measurements. 



2 Techniques 



2.1 General background 



The general principles of DFT and the pseudopotential method have been described in the 
literature p8|-|32|. Within the pseudopotential approximation only valence electrons are repre- 
sented explicitly in the calculations, the valence-core interaction being represented by non-local 
pseudopotentials which are generated by first principles calculations on isolated atoms. The 
calculations are performed using periodically repeating geometry, the occupied orbitals being 
expanded in a plane wave basis. This expansion includes all plane waves whose kinetic energy 
Ek = h 2 k 2 /2m (k the wavevector, m the electron mass) is less than the cut-off energy E cut , 
chosen so as to ensure convergence with respect to the basis set. 

In the present work the self-consistent ground state of the system was determined using a band- 
by-band conjugate gradient technique to minimize the total energy of the system with respect to 
the plane- wave coefficients. Equilibrium positions of ions were determined by a steepest descent 
method. The calculations were performed using the CETEP code (the parallel version of 
the serial code CASTEP ||32||), running on the 64-node Intel iPCS/860 machine at Daresbury 
Laboratory. 
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2.2 Generalized Gradient Corrections 



Standard DFT calculations make use of the LDA, which treats the electron density as locally 
uniform. Although this is highly successful for many purposes, it does not give accurate results 
for energy differences involving changes of bonding, which are of interest here. In the last few 
years, methods have been developed for improving the LDA through 'generalized gradient' 
approximations (GGA). It has been shown that these corrections lead to improvement in the 
calculation of total energies of atoms and molecules p7|p4| -p^1, cohesive energies [^0|-|42|, and 



the energetics of molecular adsorption and dissociation on metal surfaces We have 

recently reported a comparison of the influence of two popular GGA schemes - those due to 
Perdew and Wang P7|J5B| and Becke and Perdew [2E,27|) - on the bulk and surface properties 
of Ti02 and Sn02 
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We found that the Becke-Perdew method gives better equilibrium 
lattice parameters, while giving essentially the same surface properties as Perdew- Wang, and 
we employ Becke-Perdew in the present work. We use the technique of White and Bird [ffH] 



for incorporating GGA within the DFT-pseudopotential calculation. In some places, we shall 
also report LDA results for comparison; these were obtained using the Ceperley- Alder (CA) 
exchange-correlation function |49| . 



2.3 Generation of pseudopotentials 



First-principles, norm- conserving pseudopotentials in Kleinman-Bylander representation 
were generated using the optimization scheme of Lin et al. |5T| in order to reduce the re- 



quired value of the plane- wave cut-off E cut . The pseudopotentials used in GGA calculations 
were constructed in a consistent way by including GGA in the generation scheme. The Sn 
pseudopotential was generated using the 5s 2 5p 2 configuration for s- and p-wave components, 
and the 5s 1 5p°' 5 5<i ' 5 configuration for the d-wave. The core radii were equal to 2.1, 2.1 and 
2.5 a.u. for the s, p and d components respectively. The Ti pseudopotential was generated using 
the 4s L85 3<i 2 configuration for the s and d waves and the 4s 1 4p a5 3<i a5 configuration for the p 
wave, with core radii of 2.2, 1.5 and 2.4 a.u. for s, p and d waves respectively. The oxygen 
pseudopotential used in our LDA calculations was generated using the 2s 2 2p 4 configuration for 
the s and p waves and the 2s 2 2p 2,5 3<i ' 5 configuration for the d wave, with a single core radius 
of 1.65 a.u. For the gradient-corrected oxygen pseudopotential, we have used the single con- 
figuration 2s 2 2p 3 5 3<i ' 45 and the same core radius. The use of a core radius of 1.65 a.u. means 
that there is an appreciable overlap of the oxygen and metal core spheres in the Sn02 and 
Ti02 crystals, and in principle this could cause inaccuracies. However, direct comparisons of 



the present results with our earlier work on Sn02 ||25|| , which employed an oxygen pseudopo- 
tential with the smaller core radius of 1.25 a.u., show that any errors due to core overlap are 
very small. 

The calculations have been done using a plane wave cut-off of 600 eV for Sn02 and 1000 eV for 
Ti02- Our tests show that with these cut-offs the energy per unit cell is converged to within 
0.2 eV. 



4 



Fig. 1. The unit cell of the rutile crystal structure. Cations are represented by small black circles and 
oxygen atoms by large white circles. 

Fig. 2. Atomic structure of the clean (110) surface. 6- and 5-fold coordinated cations are represented 
by small black circles and small black circles with white centers respectively. In plane oxygen atoms 
are represented by large white circles and bridging oxygens by large gray circles. 

2-4 Densities of States 



Electronic densities of states (DOS) associated with the ground state were calculated using 
the tetrahedron method []52|,p3|, with fc-point sampling corresponding to 750 tetrahedra in 
the whole Brillouin zone. In addition, local densities of states (LDOS) were calculated by 
taking contributions only from chosen regions of real space - in practice we have used spheres 
(radius 1.5 A) centered on chosen surface oxygen atoms. For presentation purposes, we have 
broadened the calculated DOS and LDOS by Gaussians of width 0.5 eV. In the ground state 
calculations the Brillouin zone sampling is performed using the lowest order Monkhorst-Pack 
set of k points [54], as in our previous work on Sn02 |24|.|25|| . 



3 Tests on Sn0 2 , Ti0 2 and H 2 



3.1 Perfect crystals and (110) surfaces 



A more detailed discussion of the influence of gradient corrections on the calculated parameters 
of the Sn0 2 and Ti0 2 perfect crystals and (110) surfaces is given in a separate paper pETf . Here 
we present the system used in the calculations and recall the principal results obtained using 
the Becke-Perdew GGA scheme. 

The six-atom rutile unit cell of Sn0 2 and Ti0 2 is shown in Fig. 1. The equilibrium structure has 
been determined by relaxation with respect to the lattice parameters a and c and the internal 
parameter u. The equilibrium values of these parameters calculated using the Becke-Perdew 
form of GGA are given in Table 1. The agreement of a and c with experiment is very satisfactory 
for Sn0 2 and acceptable for Ti0 2 ; the values of u are excellent in all cases. 

Our calculations on the stoichiometric (110) surface of the materials have been done with the 
usual repeating slab geometry. The rutile structure can be regarded as consisting of (110) planes 
of atoms containing both metal (M) and oxygen (O) atoms, separated by planes containing 
oxygen alone, so that the sequence of planes is O - M 2 2 - O - O - M 2 2 - O etc. The entire 
crystal can then be built up of the symmetrical 3-plane O - M 2 2 - O units. The slabs we use 
contain three of these units, and our repeating cell contains 18 atoms (6 M and 12 O). The 
perfect (110) surface consists of rows of bridging oxygens lying above a metal-oxygen layer (see 
Fig. 2). The vacuum separating the slabs has been taken wide enough to ensure that interactions 
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Table 1 

Comparison of theoretical and experimental [55] values of lattice parameters a and c and the internal 
coordinate u of Sn0 2 and Ti0 2 . The theoretical values are calculated using the Becke-Perdew form of 
GGA. 



Sn0 2 Ti0 2 
calc. expt. calc. expt. 



a (A) 


4.809 


4.737 


4.747 


4.594 


c(A) 


3.159 


3.186 


3.039 


2.958 


c/a (A) 


0.657 


0.673 


0.640 


0.644 


u(A) 


0.307 


0.307 


0.305 


0.305 



Fig. 3. Calculated valence band DOS for clean (110) surfaces of Ti0 2 and Sn02- The dashed line 
represents the LDOS on the bridging oxygen site Oj. For the sake of presentation the local contribution 
has been scaled by a factor of 5. 

between neighboring slabs are small even when adsorbed H 2 and OH - are present. The width 
we use corresponds to two O - M 2 2 - O units, and is such that planes of bridging oxygens on 
the surfaces facing each other across the vacuum are separated by about 6.8 A. This vacuum 
width is somewhat greater then was used in our earlier calculations on Sn0 2 (110) p5| . (It 
will be convenient in the following to specify slab thickness and vacuum width in terms of the 
equivalent width of O - M 2 2 - O units.) 

The surface structure has been determined by relaxing the entire system to equilibrium. As in 
our previous work on Sn0 2 (110), and the work of Ramamoorthy et al. on Ti0 2 (110) ||22|| , 
we find displacements of the surface atoms of order 0.1 A, with 5-fold and 6- fold coordinated 
metal atoms moving respectively into and out of the surface, in-plane oxygens moving out 
and bridging oxygens moving very little. The modifications of the bond lengths between the 
surface atoms with respect to the perfect crystal for Becke-Perdew GGA calculations are given 
in Table 2. 

We find that the relaxed surface energy of Sn0 2 (110) is 1.16 Jm -2 and 0.84 Jm -2 for Ti0 2 (110). 
These energies are, as already disscused in Ref. ^7|, substantially lower than the corresponding 
LDA results. 

The electronic DOS and the LDOS on bridging oxygen calculated for the valence band of the 
slab systems using the Becke-Perdew form of GGA are shown in Fig. 3. The main features of 
the calculated surface DOS of Sn0 2 have already been discussed in our previous paper [pS]] . We 
note particularly the splitting off of a narrow band of states at the top of the 0(2s) band and 
the appearance of a sharp peak at the top of the 0(2p) valence band. Both these features are 
associated with bridging oxygens. The features are also present for Ti0 2 , although the splitting 
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Table 2 

Calculated bond length modifications on Sn0 2 (HO) and Ti02 (HO) with respect to the bulk values 
for Becke-Perdew GGA exchange-correlation. For atom indexing see Fig. 2 



SnOs Ti0 2 



Oi- 


Mi 


-3.8 


% 


-5.5% 


Oiv 


-Mn 


-4.2 


% 


-5.6% 


On 


-Mn 


-1.2 


% 


-1.2% 


On 


- Mi 


2.9 


% 


2.8% 


Oin 


- Mi 


4.8 


% 


4.5% 



is less marked. 



3.2 H2O and OH molecules 



We have calculated the equilibrium geometry and electronic structure of the H 2 and OH - 
molecules. Since our techniques require us to use periodic boundary conditions, the calcula- 
tions are actually performed on periodic arrays of molecules, with the repeating cell chosen 
to be a cube of side R. For H 2 0, we find that the choice R = 7 A is large enough to render 
the interactions between periodic images negligible. Matters are not so simple for the OH - 
molecule, since it carries a net charge. To make mathematical sense of calculations in which the 
repeating cell is charged, it is essential to introduce a uniform compensating background, whose 
charge density is chosen so that the net charge in the unit cell vanishes. This technique is well 
established, and has often been discussed in the literature. However, as discussed by Leslie and 
Gillan ||56|| , the total energy still converges very slowly with increasing cell size, and the leading 



term in its deviation from the asymptotic value is —a/R, where a is the Madelung constant 
for the appropriate periodic array of point charges. This leading correction can therefore be 
subtracted exactly, and we are left with a total energy which converges reasonably quickly to 
the value for the molecule in free space. In our calculations on OH - , we find that the corrected 
total energy is converged to better than 0.1 eV for R = 13 A, and the equilibrium bond length 
is converged well before this. Comparison of the calculated and experimental structural param- 
eters of H 2 and OH - (Table 3) confirms that the molecules are described accurately by the 
present methods. 

Since we are concerned with dissociation of H 2 in this paper, we have tested the influence of 
the GGA on its dissociation energy. The reaction of interest is H 2 — > OH - + H + . To compare 
our results with experimental data, we start from the experimental dissociation energy for the 
reaction H 2 — > OH + H, which is 5.11 eV. We add to this the ionization energy of H (13.60 eV) 
and we subtract the electron affinity of OH (1.83 eV), to obtain 16.88 eV. Finally, addition of the 
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Table 3 

Comparison of theoretical and experimental values of bond lengths do-H and bond angle ZH-O-H of 
H2O and OH - . The theoretical values are calculated using the Ceperley- Alder form of LDA (CA), 
and the Becke-Perdew form of GGA (BP). Experimental values are taken from Refs. |CT , £8| 







d - 


H (A) 


IH-O-H 


H 2 


CA 


0.977 


(2.1%) 


104.8 ( 0.3%) 




BP 


0.975 


(1.9%) 


103.6 (-0.9%) 




expt. 


0.957 




104.5 


OH 


CA 


0.980 


(1.0%) 






BP 


0.979 


(0.9%) 






expt. 


0.970 







zero-point vibrational energy of H2O (0.57 eV) and subtraction of the corresponding quantity 
for OH (0.23 eV) gives us the value 17.22 eV which can be compared with the calculations. 
The calculated values with LDA and GGA are 16.6 and 16.7 eV, so that even with gradient 
corrections there is a residual error of ~ 0.5 eV. 

In order to compare with the electronic DOS reported later, it is useful to note that the single- 
particle energies of the occupied molecular orbitals (MO) of H 2 can be related to experimental 
measurements. As usual, one must be cautious about comparing Komi-Sham single particle en- 
ergies with spectroscopic energies, but it is known empirically that for occupied states the 
comparison is usually justified. We therefore compare the differences of our calculated MO en- 
ergies with the corresponding differences of measured ionization energies. In the usual notation, 
the MO states of H 2 are 2ai, 16 2 , 3ai and lb\. Our calculated separations of 2a\ — 15 2 , 16 2 — 3ai 
and 3ai — l&i (experimental values from ref. |j9] in parentheses) are 11.96 eV (13.6 eV), 3.64 eV 



(3.8 eV) and 2.12 eV (2.0 eV) respectively. Our calculated separations of 2a — 3a and 3a — In 
levels in OH - are respectively 12.2 eV and 3.5 eV. 



4 Surface adsorption of H 2 



We turn now to our calculations on dissociative and molecular adsorption of H 2 on Sn0 2 and 
Ti0 2 (110). The calculations are done using the same repeating-slab geometry used for the 
bare surface, and we need to pay attention to the effects of slab thickness and vacuum width. 
Another important technical question is whether it is better to perform the calculations with 
the particles adsorbed only on one surface of each slab or on both surfaces. We refer to these 
as the one-sided and two-sided geometries. There is a strong argument for working with the 
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Fig. 4. Atomic structure of fully hydroxylated (110) surface of rutile for a) SD, b) UD and c) SM 
adsorption geometries. Atom symbols as in Fig. 2, with adsorbed oxygen atoms represented by large 
white circles and hydrogens by small white circles. 

symmetrical two-sided geometry in which the same particles are adsorbed on opposite surfaces, 
because any possible dipole moment of the repeated cell is then eliminated, and convergence 
of the adsorption energy with increasing slab thickness is likely to be improved. We have made 
tests which confirm that this is the case, and all our calculations have therefore been made 
using the two-sided geometry. 

Our tests on the effects of slab thickness and vacuum width were performed on the Sn0 2 system 
in which H + is adsorbed on top of every bridging oxygen and OH - is adsorbed at every 5-fold Sn 
site (see Fig. 4a). In these tests, every atom in the system is relaxed to its equilibrium position. 
Using the two-sided slab geometry mentioned above, we find that increase of the vacuum width 
from two to three O - Sn 2 02 - O units changes the adsorption energy by only 0.01 eV per water 
molecule, and increase of the slab thickness from three to four units gives a change of 0.03 eV 
per molecule. If the one-sided slab geometry is used, the changes are nearly ten times as great. 
The results to be presented have all been obtained with a slab thickness of three units and a 
vacuum width of two units. We stress that in all the calculations that follow the entire system 
is fully relaxed to equilibrium. 

As noted in the Introduction, we cannot be sure in advance how H 2 will prefer to adsorb. 
To study dissociative adsorption, we begin with the symmetrical full-coverage case mentioned 
above (see Fig. 4a). We then study symmetry-lowering distortions from this configuration. A 
similar strategy is followed for the case of molecular adsorption. The effect of going to lower 
coverage is than briefly examined for the dissociative case. 

4-1 Symmetrical dissociative adsorption (SD) 

The dissociative adsorption energy is obtained from the fully relaxed total energy of the slab 
system in which H + and OH - are adsorbed at both surfaces. It is calculated in the natural 
way by subtracting this energy per unit cell from the corresponding energy for the bare slab 
plus the energy of a pair of H 2 molecules; the result is, of course, divided by two, since two 
H 2 molecules (one on each surface) are adsorbed per unit cell. As usual, a positive adsorption 
energy means that the total energy decreases when the molecule is adsorbed. In the symmetrical 
dissociative (SD) case, our LDA calculations for Sn0 2 and Ti0 2 (110) yield adsorption energies 
of 1.19 and 0.91 eV per H 2 respectively. Inclusion of gradient corrections lowers the adsorption 
energies considerably to 0.48 and 0.45 eV per H 2 0. It is clear from this that gradient corrections 
have an extremely important effect on the calculated adsorption energies. A similar effect of 
lowering of the adsorption energy by inclusion of GGA has already been reported for adsorption 
on metal surfaces leading to a substantial improvement with respect to experimental 

results. 

We find that the dissociative adsorption of water causes the equilibrium surface structure to 
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Table 4 

Calculated bond length modifications on hydroxylated SnC>2 (HO) and Ti02 (HO) with respect to the 
bulk values for LDA (CA), and GGA (BP) forms of exchange-correlation. The modifications of O-H 
bond lengths are given with respect to the free OH~ molecule. For atom indexing see Fig. 4a. 



Sn0 2 Ti0 2 
CA BP CA BP 



Oi- 


Mi 


1.2 


% 


1.5 


% 


2.5% 


2.7% 


Oiv 


-M„ 


2.2 


% 


3.4 


% 


3.8% 


4.3% 


On 


-Mii 


1.1% 


1.1% 


1.4% 


0.9% 


On 


- Mi 


-1.5 


.% 


-1.4 


:% 


-1.7% 


-1.0% 


Oin 


- Mi 


-1.5 


.% 


-1.9 


1% 


-2.7% 


-3.1% 


Oi- 


Hi 


-0.8 


;% 


-o.s 


;% 


-1.6% 


-1.5% 


A 


-H„ 


-2.2 


:% 


-2.2 


:% 


-1.7% 


-2.0% 


O a 


-Mii 


-5.5 


>% 


-6.2 


:% 


-6.9% 


-7.0% 



change significantly. The relaxations found for the clean surface are greatly reduced and in 
some cases reversed. The adsorption causes 5- and 6-fold cations to move respectively out of 
and into the surface, the in-plane oxygens to move outward, and the bridging oxygens to move 
slightly outward, relative to the relaxed clean surface. Compared with free hydroxyl molecules, 
surface O-H bonds are considerably shorter, the shortening being especially pronounced for the 
OH group attached to the 5-fold surface cation. Similarly, the bond between the OH group and 
the surface cation is much shorter than the bulk O-cation bond, resembling rather the inter- 
atomic distance in the corresponding diatomic molecule. Details of inter-atomic bond length 
changes for both materials (with respect their perfect crystal values), with and without gradient 
corrections, are presented in Table 4. The changes of O-H bond lengths are given with respect 
to the free OH~ molecule. It is clear from the table that, in contrast to the adsorption energies, 
the relaxed structure of the hydroxylated surface is little affected by gradient corrections. 

Figure 5a shows both total and local DOS calculated for the hydroxylated (110) surface of 
Sn0 2 and Ti0 2 . Two oxygen sites were studied using the LDOS: the bridging oxygen atom 
Oi and the oxygen of the adsorbed OH~ group Oa- Compared with clean surfaces, two main 
modifications are apparent: 

the contribution due to the bridging oxygen atom, situated at the top of the valence 
band for clean surfaces (and giving an extra feature above the 0(2s) band) has been pushed 
towards lower energies, and hybridizes more strongly with the bulk 0(2p) band. On the 
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Fig. 5. Calculated valence band densities of states for hydroxylated (110) surfaces of Ti02 and Sn0 2 
in a) SD, b) UD c) SM adsorption geometries. The dashed line represents the LDOS on the bridging 
oxygen site Oi, and the dotted line the LDOS on the adsorbed oxygen site Oa- For the sake of 
presentation the local contributions have been scaled by a factor of 5. 

other hand, formation of the 0-H + bond gives rise to sharp bonding 3a states below the 
0(2p) band and 2a states below the 0(2s) band. 

the adsorbed hydroxyl group gives a contribution to the valence 0(2p) band due to its 
occupied lir state. It lies above the surface VBM in the case of Sn0 2 (HO) and within 
the surface VB for Ti0 2 (HO). On the other hand, bonding 3a states can be seen as a 
narrow peak below the 0(2p) band on Ti0 2 (HO), whereas on Sn0 2 (HO) they lie within 
the 0(2p) band. For both materials, the 2a peak lies within 0(2s) band. Separations of 
2a — 3a and 3a — In are respectively 10.3 eV and 6.2 eV for Sn0 2 and 11.4 eV and 5.3 eV 
for Ti0 2 . 

It is clear from the results that, especially for the adsorption energy, gradient corrections make 
a substantial difference. All the remaining calculations are performed with the Becke-Perdew 
GGA scheme only. 



4-2 Unsymmetrical dissociative adsorption (UD) 

The SD geometry described in the previous section is not the most stable one, and a small 
displacement of the adsorbed atoms from their fully symmetric positions makes the system 
relax to a configuration of lower energy. We find that this unsymmetrical dissociative (UD) 
configuration gives an adsorption energy of 1.39 eV per H 2 for Sn0 2 (110) and 1.08 eV per 
H 2 for Ti0 2 (110), so that the breaking of symmetry yields a stabilization of well over 0.5 eV. 

The nature of the new relaxed configuration is shown in Fig. 4b. The bridging (Oi) and adsorbed 
(Oa) oxygens approach each other, and the proton (Hi) attached to bridging oxygen tilts 
towards Oa, so as to form a hydrogen bond Oa~Hi of length 1.81 A on Sn0 2 (HO) and 1.80 A 
on Ti0 2 (110). The separation between Oa and Oi oxygens is 2.78 and 2.77 A for the two 
materials. These changes induce a dilation of the Oi-Hi bond (by 5.2% for Sn0 2 and 4.4% for 
Ti0 2 ) and smaller dilations of bonds between oxygens and surface cations. 

Relaxation to the unsymmetrical configuration causes noticeable changes to the DOS and LDOS 
described in the previous section. We note (see Fig. 5b) a significant broadening of the isolated 
peak due to non-bonding Itt states of adsorbed OH - . In both materials, the 3a levels are shifted 
upwards and hybridize more strongly with the valence band. Similar changes can be observed 
in the LDOS of the Oi-Hi group, with a marked upward shift of 3a into the 0(2p) band. The 
separations of the 2a — 3a and 3a — In peaks become respectively 11.7 eV and 4.6 eV for Sn0 2 
and 13.3 eV and 3.7 eV for Ti0 2 . 
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4-3 Symmetrical molecular adsorption (SM) 



We have considered the case of symmetrical molecular (SM) adsorption in which the water 
molecule bonds by its oxygen to the surface 5-fold coordinated cation and the plane of the 
molecule is the (001) plane (Fig. 4c). We find that for both materials adsorption in this geometry 
is energetically favorable, the adsorption energies being 0.78 eV per H2O for Sn02 and 0.82 eV 
per H2O for Ti02- These are considerably larger than the adsorption energies for the SD 
geometry. (Recall that we are comparing energies calculated with the GGA.) 

We find that molecular adsorption of H 2 causes only minor changes to the surface structure. 
Relative to the relaxed clean surface, the 5-fold coordinated cation moves out by 0.05 A(Sn02) 
and 0.03 A(Ti0 2 ). The bridging oxygen moves out somewhat less and the atoms of the first 
M 2 2 atomic plane move inwards. The size of the latter effect is nearly negligible for Sn0 2 (110) 
but the displacements are as large as 0.1 A on Ti0 2 (110). There is also a significant deformation 
of the adsorbed molecule: even though the O-H bond lengths remain unchanged, the angle 
between the bonds is increased by 10%. This is mainly due to modification of the water 3ai 
orbital, because of its contribution to bonding to the surface cation. The O-cation bond between 
the water molecule and the surfaces is about 10% longer than the corresponding bond in the 
bulk crystal. 

Densities of states for molecularly adsorbed water are presented in Fig. 5c. Compared with 
the DOS of the clean surface, the features attributed to the bridging oxygens (additional peak 
above 0(2s) and 0(2p) bands) remain practically unchanged. The modification due to adsorbed 
H 2 can be seen as an isolated narrow peak below the 0(2s) and 0(2p) bands, as well as a 
two-peak contribution to the 0(2p) band. These can be attributed to 2a 1; lb2, 3ai and lbi 
molecular states respectively. The energy separations between peaks are about 11.5/2.9/3.4 eV 
for Sn0 2 and 11.5/3.4/3.2 eV for Ti0 2 , which are quite close to the free molecule results 
(12.0/3.6/2.1 eV), with the biggest modification, as expected, being the downward shift and a 
small splitting (for Sn0 2 (HO)) of the 3ai peak. 

4-4 Surface dissociation 

Of the geometries we have examined, the two most stable are the UD and the SM configurations, 
with the difference of adsorption energies being ~ 0.6 eV for Sn0 2 and 0.3 eV for Ti0 2 . In 
both materials, dissociative adsorption is favored. It thus seems possible that water adsorbed 
molecularly on the surface can spontanously dissociate. Whether or not this occurs will depend 
crucially on the existence and height of the energy barrier that has to be overcome when the 
hydrogen bond is created and when the proton migrates along it. 

To investigate this problem, we have calculated the total energy of the system for a number of 
configurations along the probable dissociation path. In practice, we have chosen the reaction 
coordinate to be the horizontal (in the surface plane) separation of the 5-fold coordinated cation 
and the migrating proton. The positions of all other atoms, as well as the vertical position of 
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Fig. 6. Dependence of adsorption energy on the reaction coordinate (see text) for water dissociation 
on Sn0 2 (110) and Ti0 2 (110). 

5-fold coordinated cation and of Hi have been relaxed for each value of the reaction coordinate. 
Results on the dependence of the adsorption energy on the Mi-Ht horizontal separation for 
both materials are displayed in Fig. 6. For both materials, we find almost no energy barrier 
which would need to be overcome on the passage from SM to UD adsorption geometries. It is 
possible that there may be a very shallow minimum at the SM geometry, but this would not 
be of any practical significance. 

4-5 Adsorption in the low adsorbate density limit 

In order to get insight into the dependence of adsorption characteristics on the adsorbate density 
we have considered the half-coverage case where there is one H2O in every two surface unit cells. 
Because of limitations on computer time, we have chosen a single adsorption geometry, namely 
the UD one, which is energetically the most favorable in the high density limit. 

The practical calculations were performed with a supercell twice as big as the previous one, 
keeping however a slab thickness of three O - M2O2 - O units and a vacuum width equivalent 
to two units. The clean surface calculations were repeated for this supercell in order to ensure 
cancellation of errors between the clean and hydroxylated slabs. The positions of all slab and 
adsorbate atoms were relaxed to equilibrium. We find that the adsorption energy is greater for 
the half-coverage case, namely 1.63 eV per H 2 for Sn0 2 (HO). Compared with full coverage 
case, the atomic structure of the surface and adsorbate are only slightly modified. 



5 Discussion 

Our calculations on the energies of different relaxed configurations indicate that dissociative 
adsorption gives the largest adsorption energy and that the most stable configuration is un- 
symmetrical. The adsorption energies at full coverage are 1.39 and 1.08 eV per water molecule 
for Sn0 2 and Ti0 2 , increasing to 1.63 for Sn0 2 at half coverage. It is important to note that 
dissociative adsorption gives the strongest binding, even though we are dealing with the perfect 
surface and no defects are involved. This appears to be an entirely geometrical effect, since our 
results are very similar for Sn0 2 and Ti0 2 , in spite of their different electronic structure. In 
fact, our entire set of results both for molecular and for dissociative adsorption, and for the 
instability of the molecularly adsorbed state, is remarkably similar for the two materials. Their 
similarity may perhaps be understood by noting that they differ mainly in their unoccupied 
electronic states. 

In the TDS experiments referred to in the Introduction |T0| , [T3"|1 , the high temperature features 
at 435 K (Sn0 2 ) and 375 K (Ti0 2 ) were attributed to desorption from the dissociated state. 
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Adsorption energies can be deduced from these desorption temperatures following the analysis 
of Redhead which expresses the desorption energy rate as the product of an effective 
vibrational frequency and a Boltzman factor exp(— E^s/ZcT), where is the adsorption 
energy. Assuming a frequency of 10 13 s _1 , we find that this analysis gives adsorption energies of 
1.1 and 1.0 eV for Sn0 2 and Ti0 2 respectively. Our calculated values are thus in fair agreement 
with experiment, though they are systematically too high. It is likely that this discrepancy is 
related to the fact that our DFT calculations - even with gradient corrections - underestimate 
the energy of the reaction H 2 — > OH + H + (see sec. 3.2). The fact that the adsorption energy 
for Sn0 2 is greater than that for Ti02 is correctly given by our calculations. 

The relation of our results to the two lower temperature peaks in the TDS spectra depends 
on how one interprets these peaks. The peak at the lowest temperature appears to arise from 
desorption of molecular H 2 when more than a monolayer is present, i.e. desorption of H 2 
bound to other H 2 0. Our calculations clearly have nothing to say about this. The intermediate 
peak at 300 K (Sn0 2 ) and 275 K (Ti0 2 ) was also attributed mainly to desorption from the 
molecularly bound state. Our calculations clearly suggest that molecular H 2 directly bound 
to the (110) surface is unstable with respect to dissociation. It is relevant to note that the 
semi-empirical calculations of Bredow and Jug ]T^] on the isolated H 2 molecule on the Ti0 2 
(110) surface found a barrier separating the molecular and dissociated states, but the height of 
this barrier was only ~ 0.2 eV. There are then two possible interpretations of the experiments: 
either the water desorbed at the intermediate temperature is not bound directly to the oxide, 
but is perhaps bound to hydroxyl groups; or there is some more stable molecularly adsorbed 
state which we have not examined. Both of these seem quite plausible. The second interpre- 
tation could be probed by making a more extensive search for stable molecularly adsorbed 
configurations, perhaps using dynamical simulations, and we hope to return to this. 

Our calculations have predicted an appreciable increase of dissociative adsorption energy with 
decreasing coverage. This is consistent with experimental findings that the coverage of dissocia- 
tively adsorbed water on Sn0 2 and Ti0 2 (HO) is less than a monolayer, and can be attributed 
to the effect of repulsion between hydroxyl groups. The same effect has been found in semi- 
empirical quantum calculations on the dissociative adsorption of H 2 on Ti0 2 |p]Jl6|. 

Calculated local densities of states for the oxygen atom Oa in the adsorbed OH~ group (ge- 
ometries SD and UD) or for the adsorbed water molecule (geometry SM) strongly resemble the 
spectra of the corresponding free molecules. The principal modification concerns the molecular 
orbital directly involved in formation of the bond with the surface, 3a for OH - and 3a± for H 2 0. 
In most cases modifications consist of a downward shift (relative to the free molecule) accompa- 
nied by a small splitting, and are consistent with UPS findings on hydroxylated Ti0 2 (HO) 0. 
It is also worth noticing, that for the more stable unsymmetrical geometry the peak related to 
nonbonding In or lb\ is considerably broadened or split by the formation of a hydrogen bond 
and hybridization with the surface 0(2p) band. On the other hand adsorption of the proton 
on the bridging oxygen introduces a strong modification to its LDOS. This is clearest for the 
symmetrical dissociative adsorption case, where the additional electrostatic field of adsorbed 
H + causes a substantial downward shift of the totality of LDOS(Oi) compared to the clean 
surface. In addition a distinct peak due to bonding O-H states appears below 0(2p) band. The 
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overall structure of the LDOS on bridging oxygen for dissociative adsorption is quite similar to 
the free molecule spectrum, although for the unsymmetrical geometry broadening and splitting 
of both features in the region of the 0(2p) band become important. 



Comparison of our results with experimental UPS spectra |7|]I3|] shows quite good agreement. 
However it should be kept in mind that experimental spectra taken at low temperatures tend 
to show molecularly adsorbed water which may be in higher adsorbed layers, rather than being 
bound directly to the surface. This could explain the much closer resemblance of measured 
spectra to those of free molecules. On the other hand, the high-temperature spectra showing a 
two-peak structure which can be assigned to OH~ molecular orbitals relate well to our LDOS 
for adsorbed oxygen. A surprising feature of the experimental results is the total absence of 
signal from OH groups formed on bridging oxygens, which according to our calculations should 
also give a characteristic, two-peak structure. 



6 Conclusions 



We have studied molecular and dissociative adsorption of water on Sn0 2 (110) and Sn0 2 (HO) 
by an ab initio, density functional approach. We found that inclusion of gradient corrections to 
the LDA noticeably reduces the adsorption energies, a tendency which has been found already 
for surface adsorption on metals. For both materials, in the full coverage regime, we find that 
both molecular and dissociative adsorption are energetically favorable. However, dissociative 
adsorpton has a substantially greater adsorption energy. Investigation of a possible reaction 
path for the surface dissociation of water shows the absence of any energy barrier. We found 
an increase of adsorption energy in the half-coverage regime, which indicates the existence of 
an effective repulsive interaction between adsorbed species and suggests that the most stable 
adsorption configurations belong to the low coverage limit. Calculated valence band DOS spec- 
tra for molecular and dissociative adsorption geometries show structures similar to these of 
free H 2 and OH - molecules. Bonding to the surface introduces small shifts and broadening 
of peaks in the region of the surface valence band. 
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